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Abstract 

We have applied the Ewald-Kornfeld formulation to a tetragonal lattice of 
point dipoles, in an attempt to examine the effects of geometric anisotropy on 
the local field distribution. The various problems encountered in the computa- 
tion of the conditionally convergent summation of the near field are addressed 
and the methods of overcoming them are discussed. The results show that 
the geometric anisotropy has a significant impact on the local field distribu- 
tion. The change in the local field can lead to a generalized Clausius-Mossotti 
equation for the anisotropic case. 
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I. INTRODUCTION 



When a strong field is applied to a macroscopically heterogeneous or composite medium, 
the induced change of the medium can lead to quite interesting behavior, both in electrical 
transport and in optical response 0. These phenomena are most pronounced in the case 
of periodic composites. In the case of a cubic lattice of dipoles, elementary (dipole lattice) 
arguments show that the correction to the Lorentz cavity field (due to all the dipoles inside 
the cavity) strictly vanishes 0. 

However, when the lattice symmetry is lowered by an external means, e.g., under the 
influence of an external force/torque, the lattice is deformed, either lengthened in one di- 
rection and/or contracted in the other direction, the correction to the cavity field will not 
vanish. We call this phenomenon the geometric anisotropy ||. The lattice deformation can 
easily be realized by the electrorheological effect on a suspension of polarized particles, in 
which the particles aggregate into anisotropic structures. In this work, we will apply the 
Ewald-Kornfeld formulation HJ5[ to a tetragonal lattice of point dipoles, in an attempt to 
obtain a convergence of the infinite lattice sum and hence examine the effects of geometric 
anisotropy on the local field distribution. 

The plan of the paper is as follows. In the next section, We will consider a tetragonal 
lattice of point dipoles and apply the Ewald-Kornfeld formulation to calculate the local field. 
In section III, we perform numerical calculation and the various problems encountered in the 
computation of the conditionally convergent summation of the near field will be addressed. 
In section IV, we discuss the relation between the present formulation with the established 
local field concepts. We will show that the change in the local field can lead to a generalized 
Clausius-Mossotti equation valid for the anisotropic case. Discussion on related problems 
will be given. 
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II. FORMALISM 



Consider a tetragonal lattice with lattice constant qa along the z-axis and lattice constant 
aq~ l l 2 along the x and y axes. In this way, the degree of anisotropy is measured by how 
q is deviated from unity and the uniaxial anisotropic axis is along the z-axis. The lattice 
constants have been chosen so that the volume of the unit cell V c = a 3 remains unchanged 
as q varies. The lattice vector is given by 

R = a(q~ 1/2 l£ + q~ 1/2 my + qnz), (1) 

where /, m, n are integers. There are N point dipoles located at in a unit cell. The 
local electric field Ej at a particular point dipole at Tj can be expressed as a sum of the 
electric field of all dipoles at r Rj -: 

e, = EE t ^-p^ (2) 

j R 

where "prime" denotes a restricted summation which excludes j — i when R = and 

is the dipole interaction tensor. Eq.fl^) can be recast in the Ewald-Kornfeld form [f|.^|: 
Pi • E ; = EE [ - (P< ' Vj)B{r inj ) + (pi • roy)(Pj • r iRj )C(r iRj )] 

3 R 

- Y E ^ ex P (~^2 J ' G ) ex P( iG ' r i) E(Pi " G ) exp(-iG • r,-)] 

+ 4 M (4) 
where ray = |r« — rRj|, 77 is an adjustable parameter and G is a reciprocal lattice vector: 

27T 

G = — (g 1/2 wx + g 1/2 wy + g _1 wz). (5) 

The -B and C coefficients are given by: 

erfc(^r) 2t/ 2 2 

5 r = ; — + -7=^oexp(-r/ r , 6 

3erfc(r/r) / 4t/ 3 677 \ 2 2 

C(r) = - + + exp -77 r , 7 



where erfc(r) is the complementary error function. Thus the dipole lattice sum of Eq.(0) 
becomes a summation over the real lattice vector R as well as the reciprocal lattice vector 
G. Here we have considered an infinite lattice. For finite lattices, one must be careful 
about the effects of different boundary conditions ||. We should remark that although 
a tetragonal lattice is considered, Eq.(^j) is applicable to arbitrary Bravais lattices. The 
adjustable parameter 77 is chosen so that both the summations in the real and reciprocal 
lattices converge most rapidly. In what follows, we will limit ourselves to one dipole per unit 
cell, and the Ewald-Kornfeld summation [Eq.(Q)] can be carried out. We will consider two 
cases depending on whether the dipole moment is parallel or perpendicular to the uniaxial 
anisotropic axis. In both cases, we will compute the local field as a function of the degree 
of anisotropy q. 



III. NUMERICAL RESULTS 

Consider the longitudinal field case: p = pz, i.e., the dipole moments being along the 
uniaxial anisotropic axis. The local field E at the the lattice point R = reduces to 

E, = p g l-B(R) + n'fCW - f E II exp (^) + g, (8) 

and E x = E y = 0. The local field will be computed by summing over all integer indices, 
(I, m, n) ^ (0, 0, 0) for the summation in the real lattice and (u, v, w) ^ (0, 0, 0) for that in 
the reciprocal lattice. Because of the exponential factors, we may impose an upper limit 
to the indices, i.e., all indices ranging from — L to L, where L is a positive integer. For 
q ^ 1, the regions of summation will be rectangular rather than cubic in both the real and 
reciprocal lattices. The computation has been repeated for various degree of anisotropy with 
q ranging from 0.5 to 2.0. A plateau value for E z is found for each q within a certain range 
of r) values, indicating that convergence of the local field has indeed been achieved with the 
upper limit L — 4. 

For a larger anisotropy, however, there have been two problems associated with the 
computation of the conditionally convergent summation: (1) the range of r] that gives the 
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plateau value shrinks as q increases. Even no plateau value could be observed, say, for 
q — 0.1. In this case, the summation may still converge but possibly for a much larger L 
and the computation time may be prohibitive. (2) The local field can be used to evaluate 
the depolarization factor 0, defined by: E far = 47T0P = Ei oca i — E near , where P = p/V c is 
the total dipole moment per unit volume. The near field is the 77 — > limiting value of the 
short-range part of the summation. We find that if a direct summation over the near field 
is performed, fluctuates seriously with the increase of L, which is unacceptable. 

In order to overcome these problems, the region of summation has been taken to be 
inside a sphere of radius J a in the real lattice, and that of radius 2nJ/a in the reciprocal 
lattice, where J is a positive integer. All the contribution from the dipoles outside the 
sphere will be discarded. In this way, those dipoles that contribute significantly to the local 
field but were not considered in the rectangular box are included in the summation. As a 
result, the computation time can be shortened as a much smaller value of J can be used for 
convergence. The summation is repeated with increasing J for convergence. We find that 
the summation indeed converges to a plateau value within a wide range of r\ even for large 
anisotropy. 

The second problem is also overcome by the summation over a sphere. Although there 
are still some fluctuations of the value, the amplitude of fluctuation is greatly reduced. 
We find that converges to 0.33 at J = 8. Physically, it reminds us that is exactly equal 
to 1/3, independent of the degree of anisotropy. This implies that the far field is always 
equal to 47rP/3, as from far away, the lattice structure is irrelevant. What concerned us is 
just the total dipole moment per unit volume. The fluctuation around = 1/3 is attributed 
to the slow convergence and rapid fluctuation of the near field, rather than the local field 
which converges more rapidly. This analysis thus provides us an accurate means of finding 
the near field by subtracting the far field, i.e., 47rP/3, from the local field. 

For the transverse field case in which the dipole moments are perpendicular to the uni- 
axial anisotropic axis, Eq. (8) can still be applied to evaluate the local field by modifying G z 
to G x while taking the gradient along the direction of the dipole, say the x-axis, and obtain 
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the expression of the local field. 

The results of the local field strength (normalized to AttP/3) against log 10 g for the 
longitudinal and transverse field cases are plotted in Fig. 1(a) and Fig. 1(b) respectively. For 
comparison, the near field is also plotted on the same figure. The near field vanishes at 
q = 1, in accord with the previous result. As q decreases, the local field for the longitudinal 
field case increases rapidly while that for the transverse field case decreases rapidly. In both 
cases, when q deviates from unity, the effect of geometric anisotropy has a pronounced effect 
on the local field strengths. 

IV. CONTACT WITH MACROSCOPIC CONCEPTS 
A. Generalized Clausius-Mossotti equation 

Our present theory is of microscopic origin, in the sense that we have computed the 
lattice summation by the Ewald-Kornfeld formulation. We have not invoked any macroscopic 
concepts like the Lorentz cavity field m the calculations. However, to corroborate with 
these established concepts can lead to a modification of the Clausius-Mossotti equation valid 
for the anisotropic case. 

More precisely, we use the result of the local field to evaluate the effective polarizability 
a c s of the dipole lattice, which is given by : 

a 

where a is the polarizability of an isolated dipole, and (3 = E/P is the local field factor. 
Note that j3 = 47r/3 when q — 1. To see this, the total field acting on a dipole is the sum of 
the applied field E and the local field due to all other dipoles, hence 

p=a(E + /3P), 

where E is the applied electric field. Let P = p/V c , the above equation becomes a self- 



6 



consistent equation. Solving yields 

The effective dielectric constant e c s is given by 1 + Aira e s/V c . For a cubic lattice, (3 = 47r/3, 
e e g satisfies the well-known Clausius-Mossotti equation: 

e eff + 2 ~ 3K ' 1 j 

Thus Eq.(9) represents a generalization of the Clausius-Mossotti equation to the anisotropic 
lattice. 

The result of the effective polarizability is plotted against log 10 g in Fig. 1(c) for the 
longitudinal field case while in Fig. 1(d) for the transverse field case with a = 0.001,0.01 
and V c = 1. For a small a = 0.001, the effective polarizability is almost independent of q. 
However, for a larger a = 0.01, the effective polarizability exhibits similar behavior as the 
local field. As q decreases, a e g for the longitudinal field case increases rapidly while that for 
the transverse field case decreases rapidly. Again, when q deviates from unity, the effect of 
geometric anisotropy has a strong impact on the effective polarizability. 

B. Onsager reaction field 

A problem with the Lorentz theory is in regard to its generalization to polar media with 
permanent dipole moments fi. A simple replacement of the polarizability a by a + /i 2 /fceT 
in the Clausius-Mossotti relation leads to a divergent dielectric constant, a phenomenon 



known as the polarization catastrophy. In 1936, Onsager |10[ resolved this problem by 
introducing a reaction field: while the Lorentz local field at the origin is due to all dipole 
moments of the lattice in the absence of the dipole moment at the origin, the reaction field 
at the origin arises from the additional polarization of the surrounding dipole moments due 
to the dipole moment po at the origin. 

In the present case, the Onsager reaction field at the origin is given by: 
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R = a T r : T R0 • p , (11) 

R^O 

where a is the bare polarizability and T is the dipole interaction tensor 

Tor = VrVr-^ = T R0 . (12) 

Unlike the summation of the local field, the infinite sum [Eg. (|TTD1 for the Onsager reaction 
field is indeed absolutely convergent. No Ewald-Kornfeld formulation is needed because the 
product T:T is positive definite. For a simple cubic lattice of dipole moments, by summing 
over nearest neighbors, the result is: 

R = 12a Po /V c 2 , 

which is already close to the infinite lattice limit Rq ~ 16.8apo/V^ . It is instructive to 
extend the consideration to anisotropic lattices. The reaction field is conveniently expressed 
as Rq = Xapo/Vc, where A is the reaction field factor. Again, the summation over sphere 
has helped the convergence. We perform the direct summation in the real lattice to obtain 
the reaction field factor for various q ranging from 0.5 to 2.0. In Fig.2(a) and (b), we plot A 
versus log 10 q for the longitudinal and transverse field cases respectively. There is a minimum 
A around q = 1 (but not exactly at q = 1) in both cases and A increases rapidly when q 
deviates from unity. As evident from Fig. 2, the effect of geometry anisotropy has a strong 
impact on the reaction field strength. 

V. DISCUSSION AND CONCLUSION 

Here a few comments on the results are in order. We should remark that the present 
work employs the electrostatic (dipole) approximation. While such an approximation is 
simpler to implement and suffices in many cases, for optical properties of composites of 
metallic particles embedded in a dielectric host medium, however, one must go beyond the 
electrostatic approximation |p"T|. 



8 



It is possible to extend the present theory to statistical geometric anisotropy [fL2| , e.g., 



that being induced by the electrorheological effects . We may also extend the formulation 
to a randomly dilute lattice of dipoles to mimic mixed magnetic and nonmagnetic ions. 

Although we have employed a lattice structure in the present work, our formalism can 
readily be applied to an assembly of randomly placed dipoles by considering a more complex 
basis in a unit cell. Such an extension will be useful for the study of dielectric liquids under 
an intense electric field. 

In conclusion, we have applied the Ewald-Kornfeld formulation to a tetragonal lattice of 
point dipoles to examine the effects of geometric anisotropy on the local field distribution. 
The various problems encountered in the computation of the conditionally convergent sum- 
mation are addressed and the methods of overcoming them are discussed. We suggest that 
the large value of the derivative of the local field with respect to the degree of anisotropy 
offers potential applications as artificial piezoelectric materials as one can change the degree 
of anisotropy easily in a suspension. 
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FIGURES 

FIG. 1. (a) The normalized local field strength plotted against log 10 q for dipole moments along 
the uniaxial anisotropic axis. When q deviates from unity, the effect of geometric anisotropy has a 
strong impact on the local field strengths, (b) Similar to (a), but for dipole moments perpendicular 
to the uniaxial anisotropic axis, (c) The effective polarizability plotted against log 10 q for dipole 
moments along the uniaxial anisotropic axis. When q deviates from unity, the effect of geometric 
anisotropy has a similar effect on the effective polarizability. (d) Similar to (c), but for dipole 
moments perpendicular to the uniaxial anisotropic axis. 

FIG. 2. (a) The reaction field factor plotted against log 10 q for dipole moments along the uni- 
axial anisotropic axis. When q deviates from unity, the effect of geometric anisotropy has a strong 
impact on the reaction field strengths, (b) Similar to (a), but for dipole moments perpendicular 
to the uniaxial anisotropic axis. 



11 



Fig.l/Lo, Wan, Yu 



s 


-J 



10 



5 - 



— local field 

- near field 



-5 
-0.4 




-0.2 0.0 0.2 0.4 

log q (base 10) 



-a 



0.020 



& 0.015 



§ 0.070 

> 
■« 

u 

§ 0.005 1 



0.000 



(c) 



a = 0.01 
a = 0.001 



-0.4 -0.2 0.0 0.2 

log q (base 10) 



0.4 



-a 



8 

o 

-1 




local field 
near field 



-6 



-0.4 -0.2 0.0 0.2 

log q (base 10) 



0.4 



0.010 




0.005 



a = 0.01 
a = 0.001 



0.000 1 — 
-0.4 



-0.2 0.0 0.2 

log q (base 10) 



0.4 



600 



500 



| 400 

| 300 

s 

I 200 

a 



100 




-0.4 



Fig.2/Lo, Wan, Yu 
150 




-0.2 0.0 0.2 

log q (base 10) 



0.4 



R 



100 



X3 50 - 





-0.4 



-0.2 0.0 0.2 

log q (base 10) 



0.4 



12 



